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Abstract 

The nonlinear development of inviscid Gortler vortices in a three-dimensional boundary layer is 
considered. We do not follow the classical approach of weakly nonlinear stability problems and 
consider a mode which has just become unstable. Instead we extend the method of Blackaby, 
Dando & Hall (1993), which considered the closely related nonlinear development of distur- 
bances in stratified shear flows. The Gortler modes we consider are initially fast growing and 
we assume, following others, that boundary-layer spreading results in them evolving in a linear 
fashion until they reach a stage where their amplitudes are large enough and their growth rates 
have diminished sufficiently so that amplitude equations can be derived using weakly nonlin- 
ear and non-equilibrium critical-layer theories. From the work of Blackaby, Dando & Hall 
(1993) it is apparent, given the range of parameters for the Gortler problem, that there are 
three possible nonlinear integro-differential evolution equations for the disturbance amplitude. 
These are a cubic due to viscous effects, a cubic which corresponds to the novel mechanism 
investigated in this previous paper and a quintic. In this paper we shall concentrate on the two 
cubic integro-differential equations and in particular on the one due to the novel mechanism 
as this will be the first to affect a disturbance. It is found that the consideration of a spatial 
evolution problem as opposed to temporal (as was considered in Blackaby, Dando & Hall 1993) 
causes a number of significant changes to the evolution equations. 
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1 Introduction 


Since the work of Gortler (1940) there has been considerable interest in the boundary layer 
instability mechanism named after him. Much of the early work on Gortler vortices was shown 
to be flawed, by Hall (1982a, b), because it invoked the parallel flow approximation and thus 
ignored the effects of boundary-layer growth. Hall (1983) went on to show that for Gortler 
vortices of order-one wavenumber the ideas of unique neutral curves and growth rates are 
untenable. The stability properties of such modes depend upon the initial form and location of 
the disturbance. 

Since the early 1980’s there have been numerous theoretical studies on Gortler vortices (the 
reader is referred to the review papers by Hall 1990 and Saric 1994 for an overview of the 
subject). Most of these papers have concentrated on two-dimensional boundary layers but 
obviously many practical situations in which Gortler vortices arise will be three-dimensional in 
nature. One example of particular interest has been caused by the development of laminar flow 
control airfoils which have areas of concave curvature on the underside of the airfoil near the 
leading edge. When the wing is swept the flow becomes fully three-dimensional. Consequently 
a number of recent studies have considered the stability of Gortler vortices in three-dimensional 
boundary layers and this three-dimensionality has been shown to have an important effect upon 
the stability properties of these vortices. 

Bassom & Hall (1991) looked at the viscous and inviscid stability problems for an incom- 
pressible boundary layer flow, which could support both Gortler and crossflow vortices, over an 
infinitely long swept cylinder. They found that, for sufficiently large values of the parameter 
representing the degree of three-dimensionality of the flow, there are no Gortler vortices present 
in a boundary layer which, in the zero crossflow case, is centrifugally unstable. The inviscid 
stability problem has been extended to compressible boundary layers by Dando (1992). Similar 
results to the incompressible case are found; three-dimensionality has a stabilizing effect on 
vortices of all wavelengths except for a band of small wavelengths where the vortices are dom- 
inated by crossflow effects and are in fact of the type considered by Gregory, Stuart & Walker 
(1955). The numerical results of Dando (1992) showed that for larger Mach numbers a larger 
crossflow was needed to completely stabilize Gortler vortices over a band of wavenumbers and 
this was confirmed by the asymptotic study of Fu & Hall (1994) who considered the hypersonic 
limit. 

One of the most interesting points to emerge from the work of Bassom & Hall (1991) and 
Dando (1992) was that in the presence of a relatively weak crossflow, Gortler vortex distur- 
bances of all wavelengths are stabilized such that the inviscid modes possess some of the largest 
growth rates whilst also being neutral at certain other wavenumbers. Furthermore their govern- 
ing equation has many similarities to the Taylor-Goldstein equation which governs the linear 
stability of stratified shear flows. In fact Blackaby & Choudhari (1993) have illustrated the close 
connection between the two problems of inviscid Gortler modes in three-dimensional boundary 
layers and modes on unstable stratified shear layers, and proposed a definition of a generalized 
Richardson number (this is the parameter which characterizes the stratification of a shear flow) 
for such centrifugally-driven instabilities. The weakly nonlinear study of Bassom & Otto (1993) 
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used a classical approach to consider the stability of 0(G$) wavenumber (here G is the Gortler 
number), essentially viscous, modes in three-dimensional boundary layers. In the present paper 
we restrict our attention to the 0(1) wavenumber inviscid modes. 

It was this close connection between the two problems which encouraged the authors, in their 
desire to develop a theory describing the nonlinear evolution of these inviscid Gortler modes, to 
initially consider the nonlinear evolution of modes on stratified shear layers in Blackaby, Dando 
& Hall (1993). In order to place the latter paper and the current piece of work in context and 
understand the theory underpinning both it is necessary to review some of the recent contribu- 
tions to nonlinear critical-layer theory (see the review papers by Stewartson 1981 and Maslowe 
1986 for a more general and complete review of the theory). Over the last couple of decades 
much attention has been focussed on the nonlinear stability of non-stratihed shear layers. In 
the case where there is no vertical density variation the linear stability of the flow is usually 
governed by the familiar Rayleigh equation (to which the Taylor-Goldstein equation reduces 
for zero Richardson number). Benney & Bergeron (1969) developed the so-called equilibrium 
critical layer theory: here the mode is treated as ‘quasi-steady’ inside the critical layer as well 
as outside it. Nonlinearity affects the jump imposed across the critical layer and hence leads 
to modified results for the neutral (equilibrium) modes. Haberman (1972) extended the theory 
to include critical layers where viscosity is also significant. Some of the early studies of non- 
equilibrium critical layers include the work of Brown & Stewartson (1978), Warn & Warn (1978) 
and Hickernell (1984). The key paper by Hickernell (1984) concerned a shear layer affected by 
Coriolis (rotational) effects; here the weakly nonlinear theory leads to an integro-differential 
equation rather than the (previously) more familiar Stuart-Watson-Landau equation with its 
‘polynomial’ nonlinear terms. In fact such integro-differential equations result naturally from 
non-equilibrium nonlinear critical layer theories when the shear layer is coupled with other 
physical factors such as, for instance, Coriolis effects (eg. Hickernell 1984; Shukhman 1991); 
compressibility effects (eg. Goldstein & Leib 1989); three-dimensionality effects (eg. Gold- 
stein & Choi 1989; Wu, Lee & Cowley 1993) and buoyancy effects (eg. Churilov & Shukhman 
1988; Blackaby, Dando & Hall 1993). However the case of a ‘simple’ shear layer, not affected 
by any additional physical factors, is a special case in the sense that it does not lead to an 
integro-differential equation; instead, Goldstein & Leib (1988), found that the nonlinear evolu- 
tion of a disturbance was governed by the full unsteady nonlinear critical-layer equations. This 
difference is due to the additional physical factors, of the former cases, resulting in stronger 
singularities of the inviscid disturbance quantities at the critical level. 

At first sight, it appears that weakly nonlinear theories can only be usefully applied to 
marginally unstable flows; they rely on small growth rates and so the unstable disturbance of 
concern must be near to a neutral state. Thus it was believed that such theories are of no use 
in describing the initial evolution of ‘far-from-neutral’ unstable modes. However, a number 
of recent studies have derived integro-differential equations, using weakly nonlinear theories, 
to describe the nonlinear evolution of (general) unstable modes on a variety of shear layers 
(see the previous paragraph). These studies are based on the idea that, in actual physical 
flow situations, shear layer spreading or other external changes would result in the otherwise 
relatively unstable modes having their growth rates diminished in real terms, so that a weakly 
nonlinear critical-layer theory becomes appropriate. The work in this paper is based on the 
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assumption that boundary layer growth acts in a similar manner to shear layer spreading. This 
theory is supported by the work of Michalke (1964), Crighton & Gaster (1976) and the excellent 
comparison with experiments recently achieved by Hultgren (1992). For further discussion of 
non-equilibrium critical layer theory the reader is directed to the reviews of Cowley & Wu 
(1993) and Goldstein (1994). 

In this study we use weakly nonlinear and non-equilibrium critical-layer theories to describe 
the spatial, nonlinear development and evolution of inviscid, unstable Gortler modes in an in- 
compressible, weakly three-dimensional boundary-layer. The theory of this paper is extendable 
to compressible boundary layers and also has obvious applications to inviscid modes in a flow 
above a heated plate, similar to those considered by Hall & Morris (1992). 

Whilst different from the approach adopted in this study, there are alternate/complementary 
nonlinear theories that have been developed recently in which two or more of the flow distur- 
bances mutually interact. Such theories generally require smaller disturbance amplitudes but 
may also need the disturbances to exist in specific ‘configurations’. These other theories are 
generally referred to as wave/ wave and vortex/ wave interactions. For a discussion of wave/ wave 
interactions and resonant triads the reader is directed to the book by Craik (1985). The ideas 
behind resonant triads and non-equilibrium critical layers have been combined in works by 
Goldstein & Lee (1992) and Wu (1992) which both consider resonant triad interactions where 
the growth rates of the disturbance are controlled by nonlinear interactions inside critical layers. 
Strongly nonlinear vortex/ wave interactions were first looked at by Hall & Smith (1991) and 
their ideas were clarified and extended by Brown, Brown, Smith & Timoshin (1993) and Smith, 
Brown & Brown (1993). Recent work has in fact shown mathematical connections between 
these different nonlinear theories. Wu, Lee & Cowley (1993) in their non-equilibrium nonlinear 
critical layer study showed that the viscous limit of their amplitude equation is the same as the 
amplitude equation obtained by Brown, Brown, Smith & Timoshin (1993) in their vortex/wave 
interaction paper. 

The rest of this paper is laid out as follows. In the next section we present some back- 
ground details of the flow concerning us in this paper, namely inviscid Gortler vortices in 
three-dimensional boundary layers. In §3 the flow outside the critical layer is considered whilst 
§4 deals with the flow inside the critical layer, concentrating on deriving the non-viscous cu- 
bic amplitude equation. In §5 we consider briefly the viscous cubic amplitude equation. In 
§6 we look at some numerical solutions for the non-viscous cubic before finally drawing some 
conclusions in §7. 


2 Inviscid Gortler Vortices in Three-Dimensional 
Boundary Layers 


At this point it is helpful to recap the scalings and arguments that lead to the governing 
equation for inviscid Gortler vortices in three-dimensional boundary layers. In this paper we 
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shall consider an incompressible flow. For a more detailed derivation of this governing equation 
see Bassom & Hall (1991) and for the compressible case see Dando (1992). 

The boundary layer considered is that of a flow over the infinite cylinder y = 0, — oo < z < oo 
so that the z-axis is a generator of the cylinder and y measures the distance normal to the 
surface. The x-coordinate measures distance along the curved surface, which has variable 
curvature (1 /m)K(x/l) where to and / are length scales. The Reynolds number, i?, Gortler 
number, G, and curvature parameter, 6, are defined by 


R = 


U^l 


G = 2 R*8 

« = I 

TO 


(2.1a) 

( 2 . 16 ) 

(2.1c) 


where is a typical flow velocity in the streamwise direction and k is the kinematic viscosity 
of the fluid. 

The Reynolds number is assumed to be large, whilst 6 is sufficiently small so that as 6 — > 0 
the parameter G is fixed and of order one. The basic three-dimensional boundary layer is taken 
to be of the form 


u = U 00 [u(X,Y), R-*v(X,Y), R-i\*w(X,Y)].ll + 0(R-i)] 


( 2 , 2 ) 


with 


X = xjl Y = R 2 y/l , (2.3a, 6) 

where the parameter A* is a measure of the relative strength of the crossflow present. The basic 
state is perturbed by writing 

u = U co [u + eU(X, Y)E, vR + eR~*V(X, Y)E, 

R- 2 \*w + eR- 2 W{X,Y)E}.[l + 0{R- 2 )} , (2.4a) 

p = p{X,Y) + R- 1 P{X,Y)E , (2.46) 

where e< 1 and 

E = exp{iaR 2 z/l} . (2-5) 

We now consider the inviscid limit of the Gortler problem by introducing a scaled spatial growth 
rate, /3, and the scalings 

[U(Xo,Y), V(X 0 ,Y), W(X 0 ,Y), P(X 0 ,Y)] = [ U(Y ), 

G*V(Y), G 2 W{Y), GP{Y )] x exp jy G* /3(X)dx} , (2.6a) 

A* = Gh , (2.66) 

which were obtained independently by Timoshin (1990) and Denier, Hall & Seddougui (1991) for 
two-dimensional problems. Here X = X 0 is the local streamwise location under consideration. 
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Upon substituting these expansions into the continuity and momentum equations and letting 
G — > oo we obtain 


dV 

f3U+ — + iaW = 0 , 

(2.7a) 

(f3u + iaXw)U + u Y V = 0 , 

(2.76) 

dP 

(f3u + iaXw)V + KuU = — — , 

dY 

(2.7a) 

( fiu + iaXw)W + A w Y V = —iaP 

(2.7<f) 

From these it can easily be shown that 


d 2 V 0 (Su YY + ia\w YY ) a 2 Kuu Y 

+ a 2 r y = o 

dY 2 ( f3u + iaXw ) ( f3u + iaXw ) 2 ’ 

(2,8) 

and this is subject to the boundary conditions U(0) = 0 and V — > 0 as Y - 

oo. This is the 


equation that controls the inviscid growth of Gortler vortices in an incompressible, weakly three- 
dimensional boundary layer and, as noted earlier, it is very similar to the Taylor-Goldstein 
equation. We note that in Bassom & Hall (1991) and Dando (1992) it was chosen to scale the 
curvature, K, out of this equation, this was possible because of the local nature of the problem 
considered. At this point it is illustrative to reproduce three figures from Dando (1992) and 
consider some solutions of equation (2.8). The plots we show are of the scaled growth rate, 
jd = against the spanwise wavenumber, a, for three values of the scaled crossflow, 

A = K~ 2 A. In Figure 2.1 we have no crossflow, in Figure 2.2 we have A = 5 and for Figure 
2.3 we have taken A = 10. For further details the reader is referred to Bassom & Hall (1991) 
and Dando (1992) but these three diagrams illustrate clearly the effect that crossflow has on 
the stability of Gortler vortices. Thus in the presence of a relatively weak crossflow, Gortler 
vortices of almost all wavelengths are stabilized such that the inviscid modes possess some of 
the largest growth rates whilst also being neutral at certain other wavenumbers. 

The close relationship between equation (2.8) and the Taylor-Goldstein equation was con- 
sidered in detail by Blackaby & Choudhari (1993). In particular they proposed a generalized 
definition of the Richardson number for such vortex flows, namely 

j _ a 2 Ku c u Yc 

(/3u Yc + ia\w Yc ) 2 

where the subscript c’s denote evaluation at the critical level, Y = Y c , where (f3u -\- iaXw ) = 0. 
It is possible to calculate neutral curves for the Gortler modes (see Figure 1 in Blackaby & 
Choudhari 1993 and §6 of this paper) and then using the same definition of v as for the 
stratified shear layer problem, 

>/ = ±i(l-4.7)i , (2,10) 

we can calculate numerically the values of v on these neutral curves. It is found for the Gortler 
problem that v < — | and so from Blackaby, Dando & Hall (1993) we know that for this range 
of the parameter v there are three possible evolution equations; a viscous-cubic; a non-viscous- 
cubic (this is the equivalent of the novel non-viscous-cubic from this previous paper, it arose 


(2.9) 
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from a perturbation to the Richardson number and was denoted as the Ji-cubic) and a quintic. 
For our numerical calculations we shall consider a Falkner-Skan-Cooke profile for the base 
flow (see Bassom & Hall 1991 and Dando 1992). In order for a similarity variable to exist for 
the compressible problem (see Stewartson 1964 and Dando 1992) it is necessary to consider a 
model fluid and so for this earlier work we have taken a Prandtl number of unity. Here our 
consideration of an incompressible flow restricts us to effectively a Prandtl number of unity 
although we shall see later in §5 that this is a special case for the viscous-cubic. A similar 
definition of the Richardson number is available for the compressible Gortler problem where we 
again find that v < — | (see Dando 1993). 


3 Outside the Critical Layer 


In order to derive the desired evolution equations a study of the fundamental and other higher 
harmonics is necessary both inside and outside of the critical layer. The details of this study 
are dependent on the flow under consideration but the method is quite general and can be 
applied to other flows (the reader will note many similarities between the next two sections and 
the corresponding work in Blackaby, Dando & Hall 1993 for the stratified shear layer problem 
where some aspects are discussed in more detail). In this section we consider the flow outside 
the critical layer whilst the following two sections are devoted to that inside the critical layer. 


3.1 Scalings and Notation 

It is wise, before embarking on the analysis of the critical layer, to consider for a moment the 
various x-scales that we have in this problem. In addition to the non-dimensional boundary 
layer variable, X (defined by equation 2.3a), we have to consider: 

(i) the fast scale, X = G? X say, of the disturbance which was implicit in the scalings (2.6a) 
and was necessary for the derivation of the governing equation (2.8); 

(ii) the slower scale, X = G^/iX, over which the disturbance amplitude evolves (ji is a small 
parameter); 

(iii) the slowest scale x\ = /i~ 1 X (provided G~ * <C //), over which the mean flow varies. 

Due to the length of many of the equations in the critical-layer analysis we shall try to 
use as much abbreviated notation as possible. In particular the real functions q 0 , q 1 and q 2 will 
be defined by 

iq 0 = (f3uo + iaX 0 W 0 ) , (3.1a) 

iq i = (/3ui + ia[\iW 0 + A 0 uq]) , (3.16) 

iq 2 = (f3u 2 + ia[X 2 w 0 + X x wi + X 0 w 2 ]) . (3.1c) 

Note that the critical layer Y = Y c occurs when q 0 = 0. Also ui,u 2 , uq, w 2 , Ai and A 2 are real 
functions of x 1 and for example 

Hi = xtUoxiXo) , (3.2a) 
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(3.26) 


^2 = -xl^oxx(Xo) 

Similarly the perturbation to the curvature, Jlh, is also a real function of x 1 . From this point 
onwards dashes on mean flow quantities will denote derivatives with respect to Y . It is also 
implicit that all of the mean flow quantities in the critical layer analysis (§3.3, §4, §5 and §6) 
are the values evaluated at the critical layer. We shall not explicitly write a subscript c to 
denote this. 


3.2 The Solvability Condition 


In order to derive the evolution equation we expand the fundamental as 

Id = eVl 1] + tfiV / 2) + • • • , (3.3) 

where e is a small parameter characterizing the magnitude of the mode. Here is the 

neutral mode of the inviscid linear problem and Vp ; takes into account the X-dependence of 
the solution. We note that 

V t {1) = B±(X)V(Y) and V(Y) ~ \Y - Y c \^ +l/ as (Y - Y c ) -> 0 , (3.4) 


where the + and — sign on the B denote, respectively, above and below the critical layer. The 
second term in the expansion of Vi satisfies the equation 

LV} 2) = Q x , (3.5) 

/ o \ / o \ 

subject to the boundary conditions Vp (Y = 0) = 0 and Vp ; — > 0 as Y — s- oo, where the 
operator L is given by 


L 


d 2 ( 2 % a 2 KqUqu q \ 

,^ 2 v a % m 2 ). 


(3.6) 


and the right hand side, Q i, by 


Q i — 


q 1 ! d 2 c i 2 K 0 u 0 u 0 ^ a 2 K 0 u 0 u 0 fu 1 u\ K i\ ((jq + a 2 <q) 

- \ ¥ 0 + w + xT 1 + = 


% \9Y 2 (iq 0 ) 2 ) ( iq o y 


L o , 


% 


v; 


a) 


+ 


Uq 


d 2 a 2 K 0 u 0 u 0 ^ (u 0 + a 2 u 0 ) 


iq 0 \dY 2 (iq 0 ) 2 J 


iq 0 


dVl 


(i) 


dX 


(3.7) 


I 2 I I 2 I I 2 I 

The solution to this equation can be considered to be the sum, Hp ; = Vp PI + V{ CF , of a 

(2) ( 2) 

particular integral, Vp P j, and the complementary function, HpTp. As Y — Y c — > 0 


V$ F = B ± a {2 2\Y - Y c |^(1 + 0(\Y- Wr 1 )) + B ± b {2 l\Y - Y c |^(1 + 0(\Y - Wl” 1 )) , (3.8) 
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where a|± and are constants as yet undetermined. Note that if the Frobenius roots, | ± \ v\ 
differ by an integer then equation (3.8) is no longer appropriate (logarithms are needed). As 
such cases (v = |ra; m an integer) are isolated, we choose not to concern ourselves with them 
(and their immediate neighborhood) in this paper. 

A solvability condition for the above boundary-value problem is required. Note that: (i) 
the operator L is self-adjoint away from the critical level Y = Y c (where q 0 = 0), and (ii) the 
right-hand side of (3.5) is singular at Y = Y c . 

Following the method of Hickernell (1984) we derive the solvability condition by multiply- 
ing both sides of equation (3.5) by Vj^ and integrating over all Y, excluding the (sole) critical 
layer at Y = Y c . After integrating by parts; imposing the boundary conditions at Y = 0,oo; 
and the asymptotic forms of Vj^ and as F — ; > Y c , we find that 


£vf’e. = -(A’l-yt - lyvAHh = mb 2 A‘1 + bi^i\ 

jfl'f’Oi =2'^t(4+-*"‘vS 2 - ) ) , (3.9) 

where the barred integral signs denote the finite parts of these integrals and we have used the 
relationships B- = i~ 1 ~ 2v B + , B + = B (see §4.1). Manipulating the left-hand side of this 
equation we find that the solvability condition becomes 


f) Ft 

(Is - i~ Av h)-~; + (h ~ i~ Au h)B = 2 vB(b^l - m?}) , 

oX 


where A, / 2 , 1 3 and 1 4 are given by 


h = 


j-Yc 

fV(Y) 

Jo 


Uq 


d 2 a 2 K 0 u 0 u 0 \ (u 0 + a 2 u 0 ) 


iq 0 \dY 2 ( iq 0 ) 2 / 


iq 0 


V (Y)dY , 


(3.10) 


(3.11a — d) 


h = 


1 4 = 


iV) 

(3 = £v-m 

r 

roo 

l V{Y) 


a 2 K 0 u 0 u 0 ^ a 2 K 0 u 0 u 0 ^ u 1 + u\ ^ Ki \ ^ (q 1 + a 2 q 1 ) 


(JP_ ~ 2 * - - 

% V 5y2 ( l %) 2 J 

u 0 ( d 2 a 2 K 0 u 0 u' 0 \ (u'q + a 2 u 0 ) 


l*9o) 2 


Ur, 


u 0 


Ah 


% 


V(Y)dY 


iq 0 \dY 2 ( iq 0 ) 2 7 


iq 0 


V(Y)dY 


a 2 KqUqUq ^ a 2 K 0 u 0 u 0 ^u 1 u\ Ai\ (q 1 + a 2 q 1 ) 

u'o u 0 Kq 7 % 


<h / d 

% \dY 2 (iq o) 2 7 ( l %) 2 


V(Y)dY 


3.3 The Asymptotic Expansions as (Y — Y c ) — >• 0 

In terms of the critical layer variable rj = / u _1 (M — Y c ) we find that the asymptotes for V and 
$, as (Y — Y c ) — > 0, where the function $ is given by 

la 2 

*i = Vt- jj —rPt , (3.12) 

(2 “ v ) l % 
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(/ denotes the harmonic), for the fundamental, zeroth and second harmonics are 


Vi = e^ +v 


B 4 


\h+ v 


+ ... 


+ e/ii +u 


B 4 


% 


+ 


(1 - 2v) ( % 


q' 0 (l + 2v) 4 \q 0 


+5±i ?? |2+"4 2 2 + --- 


+ C/W 2 




+ 


$i = e^ +l/ [0] + 




B+\n\^ +v I 


% 


+ C/W 2 ^ 


(1 + 2u)q 0 2 u 0 2u 0 J \ 


vll&B ± 


'\~ v ) 


y 0 = e 2 fi~ 1+2l/ 
V 2 = e 2 /u~ 1+2l/ 
$ 2 = e 2 fi~ 1+2l/ 


+ 


( 1+ 2^ B 2 M -l + 2, 

(3 - 2v)iq 0 


+ • 


+ • 


+ • 


Un 


UQ 


(3.13a) 


(3.136) 

(3.14) 

(3.15a) 

(3.156) 


4 The Critical Layer 


The main purpose of this section is to derive a second relation between 6^ and b^} (the first 
being given by the solvability condition derived in the last section, 3.10) in order that we can 
obtain the desired nonlinear evolution equation. Denoting functions of rj by a hat we find that 
the governing equations in the critical layer are 

Ux + Vr, + Wz = -(i[Ux] , (4- 1 ) 

uqU x + ( r l u o T u i)Ux + u oV T (b^ow 0 + Aouq + \\Wq)Uz = — /U 2 [UUx + fj,UUx T VUq 
+WU Z ] ~ fi[(rju' 0 + ui)U x + (^?? 2 *9o + + (tju'q + u'^V] + k /i,~ 3 [U vv ] , (4.2) 

KqUqU + P v = — hK^KqUq + K\Uq + KqUi)IJ\ — fj, 1 [—KoU 2 ] , (4.3) 

uqWx + ( r l u o + u i)Wx + Aoo; o y -)- ( ri\ 0 w 0 + Aouq + \\Wq)Wz + Pz = 

-/j>~ 2 [UW x + ij,UWx + vw n + WWz\ - fj.[(r]u' 0 + ui)ITy + (^q 2 iql + 
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-\-(rj\ow 0 + Aiu> 0 + Aouq)! 7 ] + k/j 3 [WA ) ) ) ] 


(4.4) 


where we have retained on the right hand sides the leading order effects due to nonlinearity, 
viscosity and terms leading to a perturbation of the generalized Richardson number of Blackaby 
& Choudhari (1993), see equation (2.9). 

In the rest of this section we shall solve the governing equations for the relevant higher 
order terms of the harmonics. As suggested by the work of Blackaby, Dando & Hall (1993) we 
consider in this initial work two evolution equations, a non-viscous-cubic and a viscous-cubic. 
In this section we shall concentrate on deriving the non-viscous-cubic which is the equivalent 
of the Ji-cubic from this previous paper. In order to do this it is helpful to introduce the 
operator 


Ni, x = 


oX 


A - 'A 


(4.5) 


where / again denotes the harmonic. We have assumed that a// -3 <C 1, ie. viscous effects are 
not large enough to affect the operator, iV/ jX , at leading order. 

The solution inside the critical layer is again constructed in the form of a Fourier series; 
we expand the fundamental, zeroth and second harmonics, respectively, as follows 


Vi = e^ +u v} 1] + • • • + eV l+3i Wi (2) + • • • + t^ +v vi 3a) + • • • 



+e f il^V 1 {3b) + • 

•• + cV + • • • , 

(4.6a) 

% 

= cV 1+2i Wo (1) + • • 

• + eV"V 0 (2) + • • • , 

(4.66) 

v 2 

= cV 1+2i W 2 (1) + • • 

• + e 2 // 2l/ Vf } + • • • , 

(4.6c) 


and similarly for the P’ s, IF’s, P’s and $’s. These expansions are not necessarily completely 
ordered (depending on the relative sizes of e and n) and moreover we have only retained the 
terms necessary for deriving the desired evolution equation. The scalings follow directly from 
the outer asymptotes and/or by considering the process of harmonic generation. 


4.1 0(e/i2 +iy ) of the fundamental 


At this order the governing equations give 


PU[ 1] + Vff + iaWP = 0 , 

(4.7a) 

u 0 U^. + (qiq' 0 + iq^U^ + u' 0 V/ (1) = 0 , 

(4.75) 

o 

II 

^ c- 

< 

+ 

o 

IS 

(4.7a) 

u 0 W^ + (r]iq 0 + + Ao^o^i^ + iaP^ = 0 . 

(4.7<f) 
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From these four equations we find that 


= j(i . 

AVi-dh = 0 . 

and a solution of these two equations which matches to the outer solution is 
(1 + 2,y)zl^ fu 0 \" +u 


1 


v}'\x,v) = 


2T (}- U ) \% , 


dtB(X — t)t 2 ^ exp l— it 


% , <h 
q — + — 

Uq Uq 


= 0 


(4.8a) 


(4.86) 


(4.9) 

(4.10) 


For a detailed discussion of the solution of equations of the type of (4.8a,b) the reader is directed 
to Churilov & Shukhman (1988) and Dando (1993). Note however, that for the Gortler problem 
because of the range of u (v < — |) we do not have to consider the integral on the complex 
contour C as was required for the stratified shear flow problem. Matching with the asymptotes 
outside the critical layer yields 


B_ = i~ 1 - 2 l 'B + , ( B + = B) 


(4.11) 


The function V^\x ,q) has a single asymptotic representation in the lower-half plane (—7 r < 
arg q < 0) 

Vi 1] {X,q) = B{X)q^ v + 0{q- 1 * +v ) as \q\ -> oo . (4.12) 

Later we shall derive evolution equations for the amplitude, B(X) } but for the moment we can 
regard it as an arbitrary function that satisfies the requirement B(X) — > 0 as X — > — oo. This 
requirement is consistent with the initial condition used for our evolution equation (see §6), 
which itself is a result of insisting that the solution of the evolution equation matches to an 
‘earlier’ linear stage. 


4.2 0(e 2 /i 1+2iy ) of the second harmonic 

Equations (4.1-4) give at this order 

2/3N 2 (1) + Vff + 2iaWi 1] = 0 , (4.13a) 

U 0 U^ + 2(qiq' 0 + iqjuW + u' 0 V^ = -[fdU^ u[ 1] + + laW^U^} , (4.136) 

K 0 u 0 U^ + Pff = 0 , (4.13c) 

u 0 W^j + 2 (qiq' 0 + iq t )W^ + XqW^V^ + zaP 2 (1) = 
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+ , (4.13d) 

from which we can obtain the two equations 

h l+h-f’ = (1 - + 2(Vi ( , 1) Vi ( ; ) - VflV'S) , (4.14o) 


^, 1 - 2*2 ’ = 0 • (4.146) 

We note that the right-hand sides of these last two equations do not involve the conjugate of 
V} 1} and so and ^ 2 ^ have unique asymptotic representations as \r]\ — > 00 (in the lower half 
plane of complex rj). A solution of these equations which matches to the outside is 


y (1) = i 2t/ (l + 2n) 2 / up 


4n 0 r 2 (| - v) \ % 


f * Q*Q p OO P OO 

/ dt dt 1 / dt 2 B(X — t — ti)B(X — t — t 2 )(tit 2 )~^' 

/ 0 Jo Jo 


(ti - hfih + t 2 )^ +l, {2t + h + hyi- 1 ' expi -i(2t + h+ t 2 ) + 

( L u 0 u 0 

<H 1} = 0 , 


(4.15a) 


= 0 , (4.156) 

(again see Churilov & Shukhman 1988 and Dando 1993 for a discussion of the solution of 
equations like 4.14a where there is a non-zero right hand side). 


4.3 0(e 2 /i 1+2iy ) of the zeroth harmonic 

Here we find the governing equations are 



O 

II 

<^° 

(4.16a) 

+ X (1) ~ ~[ 

X ] U-ln + V-iUff - iaW^U^ + iaW { }}u[ 1] ] , 

(4.166) 


KotloUd + Pff = 0 , 

(4.16c) 

^o'WqVq ^ = — 

[-pu^w^ + pu { 2}w[ 1] + v^w^l + vl^wff] , 

(4.16c?) 


where we have used the notation for example to denote the complex conjugate of For 

our later working it is necessary to know Uq Vq ^ and From the above four equations 

and (4.7a-d) we find that 


rd 1 ) — —d T/dh 

u ° ~ anyuyi + 2v) dd ln ln ] ’ 

i o (1| = o , 


w" = (- , ) it , v-») v-») • 

0 y ia( 1 + 2v)iq' 0 ia( 1 + 2 u)a 2 K 0 u 0 J dr] lq lq ' 


(4.17 a) 
(4.176) 
(4.17c) 
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4.4 0(e 3 /i 2 + 3iy ) of the fundamental 


This is the order of the largest nonlinear term and for many critical-layer problems we would 
expect to obtain our jump here. However, the four governing equations are 


PU (2) + vi? + ia M 2) = 0 , 


(4.18a) 


"n/ i v + (mo + + y' 0 v} 2) = -{pti^uP} + puPuP + v 2 {1) u% + v$u£? + 


+K> (1) #fi ) - iaW^U?) + 2iaW^uP + laW^U^] , 

K 0 U 0 U[ 2) + Pff = 0 , 


(4.186) 

(4.18c) 


u 0 w[ | + {i]t% + iq^W^ + ^w' 0 vi 2) + iaP[ 2) = -[-0 2 ] W { }} + 2/3b1\ ) IT 2 (1) + PU^ ] W[ 1] 
+V 2 (1) W% + + iaW^wS + uiW^W^} , (4.18 d) 


and from these we can derive the equation 


jv. 1-Th = 0 


The solution of this which matches to the outside is 


$S 2) = 0 


(4.19) 


(4.20) 


and so we do not obtain a jump from this largest nonlinear term (this is exactly as expected from 
the corresponding work on stratified shear layers of Churilov & Shukhman 1988 and Blackaby, 
Dando & Hall 1993). 


2 

4.5 0(e/i2 +iy ) of the fundamental 


It is at this order that we first get terms on the right-hand sides of the governing equations 
due to the perturbation of the ‘Richardson number’. However, the situation is much more 
complicated than for the stratified shear flow case as instead of just perturbing the Richardson 
number we now have to perturb the quantities in our generalized Richardson number, equation 
(2.9). We also have X-derivatives of previous critical-layer terms and higher order corrections 
to the base flow values appearing on the right-hand sides. Neither of these two effects were 
present for the stratified shear layer case considered in Blackaby, Dando & Hall (1993) and they 
are a result of considering a spatial as opposed to a temporal evolution problem. These effects 
combined with the need to perturb u, w, A and K, cause most formulas from this point on to be 
considerably longer than their counterparts for the stratified shear flow problem. Specifically 
(4.1-4) give 

‘ (4.21a) 


pu[ 3a) + V^ a > + iaWp a ’ = ~Tl 


r(3a ) 


(3a) _ 


X) 

IX 
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(4.216) 

(4.21c) 


UqU^ + (r]i% + iq^U^ + u' 0 V) 3a ' > = 

~[(riu ' 0 + ui)U[x + (^ 2 *?o + + *? 2 )#i (1) + {Wo+^iWi 1] ] 

KqUqU[ ^ + Pj v ^ = — [(r]Kou 0 + K\Uo + AA)Ui)Ai , 

Uo W^ a) + (^ + iq^W^ + Ao^o^i (3a) + *«A (3a) = 

-[(qu'o + Hi)^^ + (^6 2 *9 o + + * 92 )^ / i (1) + (v^ow'q + AiuJq + Ao^)^ 1 ^ 

From these four equations we can obtain the equation for namely 




(1 + 2u) ( u” 

2 \u 0 u 0 J ) 


(4.21 d) 



(1 + 2z/) (u i ui 

2 \u 0 uo K 0/ 


V, 


(i) 


(4.22) 


and instead of solving this explicitly we find it easier to write the solution in terms of previously 
calculated functions; 


$j 3a) = 


% 


.(! + 2v )% 


2 u'o 



(h% 


. gi . <h u 0 u 0 

_/o "T" _/ ~T ~_/ _/ “r _ 


(1 + 2 z/)<? 0 2 q 0 ' 2q 


U r 


U 0 , 


(1 + 2v) fu[ Ui Ki 

O I ' ' ' T 7- 

2 \u t 


U 0 


K ( 


vi 1] + 


I u 0 u 0 3 u 0 


[2 viq'o V 2 u' Q 


UqC[q 


(1 + 2 v)t 


V 


(i) 


o/ J 


IX 


(4.23) 


3 

4.6 0(e/i2 _iy ) of the fundamental 

The governing equations give 


0 3h) + v) 3b) + iaW± 3b) = 0 , 

(4.24a) 

u 0 Uf^ + (rjiq' 0 + iq^U^ + u' 0 V 1 (36) = 0 , 

(4.246) 

K 0 u 0 u[ 3b) + = 0 , 

(4.24c) 

uo + (qiq’o + iqrW^ + Kw' 0 V^ h) + iaP} 3b) = 0 , 

(4.24d) 

and from these we obtain 


y,i + „i>r=i(i-2o4*r’ . 

(4.25a) 


(4.256) 
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These two equations have solutions 


y(3b) = $(36) = o ^ (4.2(3) 

and this implies that there is no linear contribution to the evolution equation from inside the 
critical layer, instead it all comes from outside the critical layer. Later we shall balance our 

(n) 

selected nonlinear term with this order and then derive our second relation involving b\_ ( and 
b^ 2 } by matching with the outside asymptotes. 

4.7 O (e 2 ii 2v ) of the second harmonic 


At this order equations (4.1-4) give 

+ Vff + ZiaWP = ~US > 

UqU^I + 2(77^0 + iq^U^ + To V" 2 (2) = 


(4.27a) 


2 pu[ 1] u[ 3a) + + V} 1] uil a) + v} 3a) U$ + iaW^U^ + laW^U^ 


(r]u 0 + Ui)U^ + (ri 2 iq " + 2 r]iq[ + 2iq 2 )U^ ) + (t]Uq + a()f/ 2 (1) 


KqUqU^ ^ + P 2 V — —[—K 0 U 1 ^ — [(r]K 0 u 0 + K\u 0 + K 0 ui)U 2 ^ , 

u 0 W^ + 2(qiq' 0 + iq 1 )W^ + Xqw'qV^ + 2iaP^ = 


3(2) _ 


.1 


r(bA( 1) 


77. uhLi 


(4.276) 

(4.27c) 


pui 1] w} 3a) + pu[ Za) w[ 1] + u[ 1] W^ + V^w}^ + + 2iaW[ 1) w[ Za) 


(i]u 0 + Ui)W^ + ( ri 2 iq 0 + 2 qiq l + 2z<2 2 )IL 7 2 (1) + (f]X 0 w 0 + Aiu> 0 + Aouq)!^ 


(4.27d) 


(2\ 

It is only necessary to hnd $ 2 ; to complete our working at the next order of the fundamental 
and from the above four equations we hnd that 




V\ 

. \% 


% (! + 2 ^) ( u o , u o 


Z \u 0 Uo Aq 


Uq Uq 


(1 - 2 v)iq 0 
2 a 2 K 0 u 0 


■ 2 hl , hd-hih!Ah 1) v' 1 ( A. 


% 

-^v 1 yi, vS , + 


Vn 


(1) 


24, >(..,>(.. . (I +!?i ) y(;»^») 


u 0 


a 2 K 0 UQ 


+ 


2V4 1, 4»<;*> - 2ti <1) #tA > 


(4.28) 


We again choose not to solve this explicitly but instead write the solution as 


$( 2) = vf) 


% 


*0 


Un 


1 

+ ^ I — ' 


Vi% 


q[ q Y (Uq u q \ (1 + 2v) 


(1 + 2v)q 0 2u' 0 2 uq) ' 2v\ (1 + 2z/)^ 2 + q 0 + 2% Vug + u 0 . 


u'n 
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u-i Ah 
+ — + T~ 

Uq Kq 


I_y(l) 1 / UqUq ou o U 0 q 0 \ ,y(l)y(l) 

[4mq 0 \ 2u 0 2 (1 + 2l/)q’ 0 J \ " [;■ I I J -'v,. ' ’ V ' J 


+V+V+ 


1 ^ UoU 0 ! 3u 0 u 0 g 0 

2?o 


1 


h+++) 


1 


+ 


L(i - 4z/ )*?o u 1 + 2z/ k 

/ uplto (3 - 2v)u 0 


0 U oJ J 


u, 


+ y (1) y (1) - 

^ lj ? 1))X 


Uq([q 


|_4z/l(l — 2i')(iq' 0 ) 2 \ 2u' 0 2(1+2;/) (1 + 2 v)q c 


u n 


_Z/(1 + 2z/)(z 9o ) 2 _ 


o/ J L 


vl 1] v {1) . - V yi ~’Vi 

1 1 r]T]X IX lr m 


Ki)t >(i) 


(4.29) 


4.8 0(e 2 /i 2iy ) of the zeroth harmonic 


The governing equations give at this order 


u 0 U^ + u 0 V 0 {2) = ~ 


T/( 2 ) — n 

v 0'0 — U 5 

JL ^^OyW + c/f^yW + u^vi 3a) + &{ 1) vlf a) 'j + 2 


(4.30a) 


+u[ 1) u^ Jt 


(f]u 0 + + (rju'o + a()y 0 (1) 


(4.306) 


KqUqUq ^ + Pq^ — — [KqU[ ^Pii] — [(tjKoUq + KqU\ + K\Uq)Uq , (4.30c) 


uoW" + Aou, 0 V 0 ( ^ — — 


d_ 

dr) 


— ( + fFT^m + w\ L ’v yoa ’ + + ( u yi, w 


r (3a)-rV(l) 


h 1 )fV( 3a ) 


dl)f/( 3a ) 


(Tw(l) 

-IX 


+u (1 )w (1 ~ } + w (1 )u (1) + ty, ( 1 ) p (1) - 

^ -! 1X^ - 1 1X^ 1 -IX 


(^^^^^^^(TyAo^o + AirCo + Ao^i)^ 1 ) . (4.30c?) 


For our working in the next subsection, at order (e 3 /r 2 + 3i/ ) of the fundamental, we need to be 
able to substitute in for 

f Q 2 « KqUq A y(2) ■ ^(2) 

^ (1-2 v)iqj ° v + ° v ’ 

and from the above four equations we find that 


(4.31) 


/3 - ^ 0 ) U$ + iaWff 


(1 - 2 u)iq 0/ 

+ (b- 


1 


d 2 


'% (1 + 2 v)uq (1 — 2v) 


Ur, 


(1 + 2 v)u 0 ) \% (l + 2i')u 0 q' 0 


v(Y + 2u)iq 0 dr] 2 [ 

'U-i Ah 

{ — + ~rr ) 
\u 0 K 0 


V\ ~ 
. \% 


2u’ 0 


2ur 


(1 + 6 u)u' 0 \ t (q[ (1 + 6 v)u' oqi (1-2 v)(u li K 1 \ (1 + 2v)u 1 A K (1) ^ (1) 


2u’ n 


— 1?7 v lr) 


. (4.32) 
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4.9 0(e 3 /i 2 + 3iy ) of the fundamental 


It is at this order that the non-viscous cubic jump will arise. The governing equations (4.1-4) 
eventually yield at this order 


where in particular 


h.i+.vf 1 = { Xp-*x\ i] + 
y = d 4) , 


4 4) = f} 2) + 4 4) 


and F[ 2 ^ is a function proportional to whilst R ^ is given explicitly by 


( 4 ) = (1 ~ 2v)(rjiq 0 + iq t ) 
3 (1 + 2v)a 2 K 0 UQ 




(1 - 2 v)ig 0 
2a 2 K 0 UQ 


(4.33a) 

(4.336) 

(4.34) 


(rjiq'o+iqjl-VffV^ 


+ Vf'Vl tm 


(1),-,(1) \ _ (3 - -2i/)i4 0 ■ (1) - (H (1 - 2i4‘<fo .-,(i),-,|i) 

g v -1 v 2r] "V 2 ^2 v -I t] 


u n u-\ Ah 

+ ( V~ T n h TT 
a 0 a 0 A 0 


- by h',!,! + vf’y 1 .!,, - ijbio;!, - A’ hh 4) + mm-?’ 


+ 




r ( 1 )<r>( 3a ) 


-*grW + iaW^) - «5>#g> - b’-'TS 1 , + T ’{ (,8 - { i“_2^Q tr S' + taW < 


(2) 
0 7 ] 


(1 - 2 v)iq 0 
2a 2 K 0 UQ 


Ai) ( ,y(i),y(i) 1 

l *lr/ *— IrjTj - 


(■-thy) + + tpA-S 


(4.35) 


To derive an evolution equation the asymptotic representation of V is needed. This is 
V} 4) = c { ±\* +v + + 0{j]-^ +v ) as T) -> Too, 

and in particular is given by 


(4.36) 


D (4) - dx = - 


( 4 ) _ 


i 2 v 


2z/u 0 T(| — v) \u 0 




[ dr] [ dt F v R^\X — t, rj) exp {— it 

•y—oo JO l 


% , ?i 

•q— + — 

Uq UqJ 


(4.37) 


After splitting R^ into three parts: firstly one containing X derivatives; secondly one involving 
terms of the second harmonic; and finally one that contains neither, we denote the pieces of 
the jump due to these parts by by subscript 1, 2 and 3 respectively so 


(Tf - D {4) ) = (D (4) - D { % + (D (4) - D { % + (D (4) - D { % , 


(4.38) 
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where the terms on the right-hand side are given in Appendix A. Matching the inner asymptote, 
(4.36), with the outer asymptote, (3.13a), leads to the relations 


flW = b<j»B + 


(4.39a, b) 

which can be combined to give 



( D S_ 4) - D ( * y ) 

II 

j33 

o-* 

1 

1 

(4.40) 

and so we have found our second relation between 6^+ and b^}. 



5 The Viscous Cubic 


In order to derive a viscous-cubic evolution equation it would be necessary to consider some 
additional critical layer terms and instead of the expansions (4.6) we would have to consider 

Vi = + • • • + eV l+3i Vi (2) + • • • + tii* +v v} 3a) + • • • + t^- v v} 3i) + • • • 

+e 3 / u“ 3 2 +3i/ I> 1 (2) + b £/«//“ l +i/ I> 1 (5) + heV^ +3 T 1 (6) + -" , (5.1a) 

V 0 = e V 1+2i % (1) + ' ' ' + eV 2i % (2) + ' ' ' + eV” 4+2i/ I>o (3) + ' ' ' , (5-16) 

V 2 = e 2 f/,~ 1+2l 'V 2 ^ 1 ' > + • • • + e 2 / u 2i/ y 2 (2) + • • • + e 2 a / u _4+2i/ y 2 (3) + • • • , (5.1c) 

We would obtain the viscous-cubic by balancing the nonlinear term at 0(e 3 « / « _ “ +3i/ ) with the 

3 

term at and then matching with the outside asymptotes in order to get our second 

relationship between 6^+ and b^}. 

However, for our current problem this proves to be impossible. It is illustrative to consider 


the governing equations in the critical layer at 0(e/c// 2 +i/ ) where viscous terms first appear on 
the right-hand sides 


pu[ 5) + Vff + iaW{ 5) = 0 , 

(5.2a) 

+ (mo + + y' 0 v} 5) = trill , 

(5.26) 

K 0 u 0 ui 5) + Pff = 0 , 

(5.2c) 

u 0 W^ + (r]iq 0 + iq^W^ + AoW o y i (5) + iaP [ 5) = W^H . 

(5. 2d) 

From these we obtain the equation 


A.i-dS 51 = 0 ■ 

(5.3) 

which has a solution 

$i 5) = 0 . 

(5.4) 
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This result is actually implied by the corresponding work on stratified shear layers (see Blackaby, 
Dando & Hall 1993 for details) because there contains the factor ( Pr — 1), where Pr is the 
Prandtl number. In this paper we have considered an incompressible fluid and hence we have 
effectively a Prandtl number of unity and as a consequence = 0. In fact the expression 
for the whole viscous cubic jump (see equations 4.45a,b and 4.46 in Blackaby, Dando & Hall 
1993) contains the factor ( Pr — 1) and hence this viscous cubic jump does not occur for our 
present incompressible Gortler problem. Obviously this nonlinear term will be important if 
we take a more realistic value (eg. 0.72 for air) of the Prandtl number but, we note that 
for the corresponding compressible problem it is necessary to consider a model fluid in order 
for a similarity variable to exist for the base flow and this nonlinear term may prove hard to 
investigate (see Dando 1995 for a discussion of this nonlinear term for the stratified shear flow 
problem). We note that for the current problem viscosity still plays a role in the sense that the 
non-equilibrium critical-layer analysis of §4 is only valid when K,fi~ 3 <C 1. 


6 The Evolution Equations 


As the viscous-cubic is not present for our problem (because the Prandtl number is effectively 
equal to unity) we can write a composite evolution equation in the form 


(h 


r'-ofl + i/, 


i- 4 v I 2 )B = 


2//H| e V 3+4i/ (^ 2 + - + cV 7+6i/ (^ 2 + - (6.1) 

where we have derived (&)+ — i Av b\_)j 1 in §4 and (&)+ — i Av b\_) q denotes the total part of the 
jump due to the quintic nonlinearity. We could derive (6^ — i~ 4 l/ b^}) q explicitly in a similar 
manner to that in which we calculated the corresponding term in our study on stratified shear 
layers. However, this would be an extremely long and complicated task for the present problem 
and for the moment we shall concentrate instead on obtaining some numerical results for the 
non-viscous cubic evolution equation (this will determine how relevant the quintic nonlinearity 
is). This evolution equation is valid provided fj, , otherwise the linear term will require 

modification (see for instance the amplitude equations 4.1a,b in Hall & Smith 1984). 

In Figure 6.1 we show the regions of validity of the various amplitude equations. The 
non-viscous cubic becomes important when the amplitude of the disturbance, A ~ H 2 2i/ (ie. 
when the cubic term becomes as large, 0(1), as the linear terms in the evolution equation, 
6.1). The quintic term is of the same size as the cubic term when A ~ fj 2 ~ v and so, if the 
cubic evolution equation results in a significant increase in the amplitude of the disturbance, 
the quintic nonlinear term will become the dominant term in the evolution equation. 
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6.1 Numerical results for the non— viscous cubic 


We now concentrate on obtaining some numerical results for the non-viscous-cubic amplitude 
equation that we derived in §4. To ease numerical calculations the jump expression is turned 
into kernel form, 


(t* 4 ) - D ( y) = 


( 4 ) _ i 1 2v (l + 2v) 


2 /tt 


2 3 z/ 2 UqF 4 (| — v) 



l-4i/ 


ds s 


2-Au 


daB(X — s)B(X — as) 


B(X — (1 a)s) + Gfi 2 (cr) 


./ \ —4b 


+ 


B(X — as)B(X — (1 + < 7 ) 5 ) G < 2 i( cr ) T G < 22 ( cr ) 


i~ 2l/ (l + 2z/) 2 q 0 

2 3 n 2 n^r 4 (| — u) \u 0 ) Jo 

i~ 2v ( 1 + 2 vf 


+ ■ 


dss 1 -^ daB(X-s) 




fCO 1 

2 — 4i/ 


ds s 2 


da 


2 8 u 2 Mq( 1 — 2n)T 2 (| — u)T(l — 2n) yn 0 
B(X - s)Bj t (X - as)B(X - (1 + <t)s)G 3 i(<t) 

+Bx(X - s)B(X - as)B{X - (1 + <t)s)G 32 (<t) 


( 6 , 2 ) 


where the kernels G'n(cr), G\ 2 (a) } G 2 i(a) } G < 22 (c r ), G < 3 i(c r ), G < 32 (c r ) are all defined in Appendix 
B. An inspection of this kernel form shows that things are considerably more complicated than 
for the corresponding temporal stratified shear flow case. The jump now has real and imaginary 
parts and there are cubic terms in the evolution equation which contain X-derivatives. The 
inclusion of spatial derivatives of the amplitude in the nonlinear terms of the evolution equation 
is a relatively novel feature for non-equilibrium critical-layer studies. Previously streamwise 
derivatives have only been seen in the study of Churilov & Shukhman (1994) on the spatial 
evolution of helical disturbances to an axial jet and spanwise derivatives in the work of Gajjar 
(1995) on stationary crossflow vortices in fully three-dimensional boundary layers. As with 
these other works the inclusion of these nonlinear spatial derivative terms is wholly associated 
with the spatial formulation of the problem (note that there were no equivalent terms in the 
related temporal stability study for stratified shear layers). It is no longer possible to determine 
whether the ‘kernel is positive or negative’ and so we are unable to deduce the behavior of the 
disturbance amplitude, as was done by Churilov & Shukhman (1988) and Blackaby, Dando & 
Hall (1993). Instead we are forced to consider numerical solutions. 

As mentioned earlier, in §2, it is possible to calculate neutral curves for the Gortler modes. 
In this paper we shall concentrate on obtaining numerical results from the amplitude equation 
for the first Gortler mode, as it becomes progressively harder to obtain neutral curves and 
n-values for the higher modes. In Figure 6.2 we have plotted the neutral curve for this first 
mode. We note that the flow is stable to this mode inside the curve plotted and that the left- 
hand branch tends towards a constant wavenumber (a = 1.305) as the value of the crossflow 
parameter tends to infinity. This limit corresponds to the neutral crossflow vortices that are 
associated with the inflection point of the directional profile, and were first investigated by 
Gregory, Stuart & Walker (1955). In Figure 6.3 we plot the values of z/, as determined from 
(2.9) and (2.10), on the left-hand branch of this neutral curve for the first Gortler mode. We 
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concentrate our attention on this area of the neutral curve as the numerical computations 
become more complicated (in particular calculating the integrals A, / 2 , I 3 and / 4 ; see equations 
3.11a-d) as v becomes increasingly negative. 

Substituting (6.2) into equation (4.41) and then matching the two relationships we have 

/ o \ / o \ 

derived between b \_ ( and b\_ (equations 3.10 and 4.41) leads to the evolution equation 


dB f f c 

at = 7lB + 72 A i 


ds s 


2-Au 


daB(X - s)B(X - (js)B(X - (1 + a)s) G n (<r) + G 12 (a) 


+74 / ds s 


l-4i/ 


/O JO 

roo rl 

+75 / ds s 2 ~ 4lJ / da 
Jo Jo 


daB(X - s)B(X - as)B(X - (1 + a)s) 


G 2 l((T ) + hr22 (^ ) 


B(X - s)B jt (X - as)B(X - (1 + (t)s)G 31 ((t) 


+B jt (X - s)B(X - as)B(X - (1 + <t)s)G 32 (<t) 
where the constants 71, 72, 73, 74 and 75 are given by 


7i = 


~ h) 


72 = 


(/ 3 - t-^h) ’ 

i~ 2v ( 1 + 2 vf 


\ —iv 


\u 0/ 


73 = 


% 


74 = FTF7 ’ 

(1 + 2 v) 

75 “ 32(1 - 2n)r(l - 2n) ' 

For our numerical calculations it is convenient to introduce the variable 


(6,3) 

(6.4a) 

(6.46) 

(6.4c) 

(6.4d) 

(6.4e) 


x = exp{27i r X} , 


(6.5) 


with 7i r denoting the real part of 71 (this is similar to the ‘logarithmic time’ variable which 
was introduced by Churilov & Shukhman 1988). We also set 

B(X) = b(X)exp{ llr X} , (6.6) 

and assume that b(X) — > 1 as X — > 00, note that the requirement B(X) — > 0 as X — > —00 is 
satisfied as 71,, is negative. After this we find that the evolution equation reduces to 

= f 1 daGUa) / 00 drr 2 - 4i/ e- 27l ’' r 6(7e- r/(1+ff) )6(7e- OT/(1+ff) W7e- r ) 

Ox Jo Jo 
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Jo Jo 

r 1 /-oo ,, N f)h( re~ ar / ( - 1+a ))_ 

-\-x / daG^A^) / c/rr 2-4i/ e -27lrr &(;re -r/ ^ 1+<7 - ) ) — b(xe ~ r ) 

Jo Jo ox 

r 1 /-00 <Wrc _, 7( 1 +Hi 

+5; / daGzBicr) / drr 2 ~ Au e~ 2llrT — &(;^e -<7r/ ^ 1+<7 - ) )&(;^e -, ’) , (6-7) 

Jo Jo ox 

with the new kernel functions being defined as 

77i = 7) 7T~ \2-4i/ [73(^11 + G12) + 7175(7731 + 7732)] , (6.8a) 

27 i r (l + a) 2 4i/ 



( 6 . 86 ) 

(6.8c) 

(6.8d) 


Finally we present some numerical solutions for the evolution equation (6.7). We note 
that we have not eliminated the parameter x 1 from the equation (see 3. 1 lb,d and 3.2) as 
it does not occur in all terms (cf. Blackaby, Dando & Hall 1993 where the parameter J\ was 
eliminated from the final evolution equation) and so for these numerical solutions we have taken 
the representative value x 1 = 1.0. We consider the effect of a small change in the crossflow 
so we have also taken Ai = 1.0 and Ki = 0. In Figures 6.4a-d we show the evolution of the 
amplitude b{x) for the case when A 0 = 10.0 (ie. close to the neutral point a = 1.371). We have 
plotted the real part of b(x) on the x-axis and the imaginary part on the y-axis and from these 
Figures we can see that the disturbance evolves in a spiral (in the complex plane of b) with 
rapidly increasing amplitude. In fact the amplitude increases so rapidly that the structure of 
the spiral is somewhat lost and this is why we show four plots of the evolution (with different 
scales). In Figures 6.5a-b we consider the evolution of b{x) when A 0 = 25.0, this is close to the 
neutral point a = 1.316. In Figure 6.5a we have plotted the argument of b(x) (in the range — 7r 
to +7r) and in Figure 6.5b we show the amplitude. The evolution again takes the form of a 
spiral in the complex plane of b{x) with a large increase in amplitude but this does not occur as 
quickly as for the previous case considered. We note that the rapid increase in the disturbance 
amplitude for these two cases means that the disturbance will enter region Illb in Figure 6.1 
where the dominant nonlinear term in the evolution equation is the quintic. This means that it 
may be necessary in the future to derive the quintic nonlinear term via a critical layer analysis 
similar to that of §4 and obtain numerical solutions for this new evolution equation. However, 
if we increase the crossflow further we find that this rapid increase in the disturbance amplitude 
does not occur. In Figures 6.6a-b we plot the argument and amplitude of b(x) for A 0 = 55.0 and 
we can clearly see that the amplitude increase has been greatly slowed. If we further increase 
the crossflow to a value of A 0 = 75.0 in Figures 6.7a-b we can see that there is hardly any 
increase in the disturbance amplitude. Consequently the disturbance will not enter region Illb 
in Figure 6.1 where the evolution equation would be dominated by the quintic nonlinearity. 
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We have shown that, at least for the nonlinear evolution equation with cubic nonlinearity 
due to the novel mechanism, crossflow has the same effect as with a linear stability study: it 
has a stabilizing influence on inviscid Gortler vortices. As mentioned earlier we have looked at 
a flow with a Prandtl number of unity, this was done because we wish to extend this work to 
the compressible regime where it is necessary to consider a model fluid in order for a similarity 
variable to exist for the base flow. However, this particular choice of Prandtl number turns 
out to be a mathematical anomally as the viscous-cubic jump expression contains the factor 
(rj — 1) and so is not existent for this problem. However, if we considered a more general 
Prandtl number then the viscous-cubic would be important because, if we consider Figure 6.8, 
which shows the regions of validity of the various amplitude equations when the viscous cubic is 
present, we see that our above results suggest that for larger values of the crossflow parameter 
the disturbance will pass from region Ilia to region IIIc where the dominant nonlinear term is 
the viscous cubic. Furthermore, for the related stratified shear flow problem (Dando 1995) it 
is known that the viscous-cubic jump results in an explosive burst-like growth. So it would 
appear to be necessary next to consider the viscous-cubic, with r] = 0.72, for the Gortler 
problem and determine whether a similar explosive growth occurs (so a disturbance with a 
reasonable crossflow moves from region Ilia to region IIIc and then onto Illb) or if sufficient 
crossflow can prevent this and hence for large values of A 0 (which still corresponds to a small 
amount of three-dimensionality) inviscid Gortler vortices will be nonlinearly stable. 


7 Conclusion 


In this paper we have considered the nonlinear development of inviscid, incompressible Gortler 
vortices in a three-dimensional boundary layer using non-equilibrium critical-layer theory. 
From the earlier work of Blackaby, Dando & Hall (1993) we knew that it was necessary to 
consider three different amplitude equations, a viscous cubic, a non-viscous cubic and a quintic. 
In this study we have concentrated our attention on the two cubic evolution equations (the non- 
viscous cubic in particular) as these will be the first to affect the disturbance amplitude. In 
§5 we showed that for the problem we have considered, with a Prandtl number of unity (a 
necessary choice in order for a similarity variable to exist for a base flow which is compressible; 
we wish to extend this work to the compressible problem), the viscous-cubic jump is zero. 
Even if this were not the case we note that, based on our assumptions concerning viscous- 
spreading effects resulting in an unstable linear disturbance mode approaching a later neutral 
state, initially our disturbance will lie in the bottom right-hand corner of Figure 6.1 or 6.8 and 
this indicates that the base evolution equation with the non-viscous cubic nonlinearity deserves 
the first attention. Consequently we have concentrated in this paper on obtaining numerical 
solutions of the non-viscous cubic evolution equation. 

Our numerical solutions show that for small values of the crossflow parameter the dis- 
turbance amplitude evolves by describing a spiral, in the complex plane, of rapidly increasing 
amplitude. This large increase in the amplitude of the disturbance means that the evolution 
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process will soon move on to a stage where the evolution of the mode is governed by an integro- 
differential amplitude equation with a quintic nonlinearity. In the related study of Blackaby, 
Dando & Hall (1993) the integro-differential equation with a quintic nonlinearity led to contin- 
ued growth of the disturbance amplitude which resulted in the effects of nonlinearity spreading 
to outside the critical level, by which time the flow has become fully nonlinear. Although 
we have pointed out the many similarities between the problems of the nonlinear evolution of 
modes on unstable stratified shear layers and the nonlinear evolution of inviscid Gortler vor- 
tices in three-dimensional boundary layers we are unable to deduce the behavior of the quintic 
integro-differential equation for the Gortler problem from the earlier results of Blackaby, Dando 
& Hall (1993) because of the major differences in the kernels for the two problems. A large 
number of these differences (including the novel feature of streamwise derivatives in the nonlin- 
ear terms) are due to the spatial formulation of the Gortler problem as opposed to the temporal 
formulation of the problem considered in Blackaby, Dando & Hall (1993). Consequently the 
daunting task of obtaining numerical solutions for the quintic integro-differential equation for 
the Gortler problem will need to be undertaken in order to clarify further the nonlinear stability 
of inviscid Gortler vortices in flows that are only very slightly three-dimensional. 

However, for larger values of the crossflow we find that the disturbance amplitude still 
evolves by describing a spiral in the complex plane but its amplitude increases only very slowly. 
Consequently when the Prandtl number is not unity the disturbance will pass from region 
Ilia in Figure 6.8 to region IIIc (as opposed to Illb when the crossflow is small) and will be 
controlled by an integro differential equation with a cubic nonlinearity due to viscous effects. In 
an investigation of the corresponding stratified shear flow problem (Dando 1995) it was found 
that the viscous-cubic resulted in large amplitude growth (of either a burst-like nature or by 
oscillating). So it would be desirable to investigate the viscous cubic amplitude equation for 
the Gortler problem (with a Prandtl number not equal to unity) and determine whether the 
same thing happens or if this amplitude growth is prevented by a sufficiently large crossflow (as 
we have found for the novel cubic amplitude equation in this paper). This may however, prove 
difficult because for the compressible problem as we need to consider a model fluid (Prandtl 
number of unity) in order for a similarity variable to exist for the base flow. 

As we have already mentioned we wish to extend the work in this paper to a compress- 
ible flow. Linear stability studies have showed that a larger crossflow is needed to stabilize 
the Gortler vortices when the flow is compressible and it would be interesting to see what ef- 
fect compressibility has on the results presented in this paper. We also note that the theory 
developed here can be applied to a nonlinear study of the inviscid vortex instabilities in the 
three-dimensional flow above a heated plate. 
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Appendix A 


The three pieces of the non-viscous cubic jump (see equation 4.38) are given by 
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where H denotes the Heaviside function. 


Appendix B 


For our numerical calculations it is necessary to convert the jump expression into kernel form, 
(6.2). The kernels in this equation are then given by 
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where F is the hypergeometric function of one variable and +1 the hypergeometric function of 
two variables (see Erdelyi 1953 and Abramowitz & Stegun 1964 for details). 
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Figure 2.1. Real parts of the scaled growth rate, /3, when there is no crossflow. 


Figure 2.2. Real parts of the scaled growth rate, /3, when the scaled crossflow, A has a value 
of 5. 
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Figure 2.3. Real parts of the scaled growth rate, /3, when the scaled crossflow, A has a value 
of 10. 
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Figure 6.1. A diagram of the various regimes of the critical layer when r] = 1. 

I: viscous, steady critical layer; Landau-Stuart-Watson equation. 

II: strongly nonlinear, equilibrium critical layer; Benney & Bergeron theory. 

Ilia: unsteady critical layer; largest term in (6.1) is cubic and due to non-viscous effects. 
Illb: unsteady critical layer; largest term in (6.1) is quintic. 

Solid lines represent boundaries with areas where the critical layers are not unsteady, dashed 
lines represent boundaries between different base evolution equations from (6.1) and the dotted 
lines indicate the threshold values at which the nonlinear evolution equations become valid. 
The thick line on the diagram indicates the evolutionary path of the disturbance for small 
values of the crossflow parameter. 
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Figure 6.2. The neutral curve of the first Gortler mode. 


Figure 6.3. The values of v for the left-hand branch of the neutral curve of the first Gortler 
mode. 
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Figure 6.4a. The real and imaginary parts of b(x), showing its evolution for A 0 = 10, part a. 


Figure 6.4b. The real and imaginary parts of b(x), showing its evolution for A 0 = 10, part b. 
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Figure 6.4c. The real and imaginary parts of b(x), showing its evolution for A 0 = 10, part c. 


Figure 6.4d. The real and imaginary parts of b(x), showing its evolution for A 0 = 10, part d. 
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Figure 6.5a. The argument of b(x) (in the range — 7r to +7r) for A 0 = 25. 
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Figure 6.5b. The amplitude of b(x) for A 0 = 25. 
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Figure 6.6a. The argument of b(x) (in the range — 7r to +7r) for A 0 = 55. 
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Figure 6.6b. The amplitude of b(x) for A 0 = 55. 
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Figure 6.7a. The argument of b(x) (in the range — 7r to +7r) for A 0 = 75. 
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Figure 6.7b. The amplitude of b(x) for A 0 = 75. 
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Figure 6.8. A diagram of the various regimes of the critical layer when r] ^ 1. 

I: viscous, steady critical layer; Landau-Stuart-Watson equation. 

II: strongly nonlinear, equilibrium critical layer; Benney & Bergeron theory. 

Ilia: unsteady critical layer; largest term in (6.1) is cubic and due to non-viscous effects. 
Illb: unsteady critical layer; largest term in (6.1) is quintic. 

IIIc: unsteady critical layer; largest term in (6.1) would be cubic and due to viscosity. Solid 
lines represent boundaries with areas where the critical layers are not unsteady, dashed lines 
represent boundaries between different base evolution equations from (6.1) and the dotted lines 
indicate the threshold values at which the nonlinear evolution equations become valid. The 
thick line on the diagram indicates the evolutionary path of the disturbance for sufficiently 
large values of the crossflow parameter. 
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